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Abstract 

We present for numerical use the analytic continuations to complex arguments of those basic 
Mellin transforms, which build the harmonic sums contributing to the 3-loop anomalous 
dimensions. Eight new basic functions contribute in addition to the analytic continua¬ 
tions for the 2-loop massless Wilson coefficients calculated previously. The representations 
derived have a relative accuracy of better than 10“' in the range x e [10“ 6 , 0.98]. 


1 Introduction 


A very efficient way to calculate the evolution of the parton densities in deeply inelastic scattering 
consist in solving the evolution equations in Mellin-iV space [1,2]. With the advent of the 3-loop 
anomalous dimensions [3,4] the scaling violations of the deep-inelastic structure functions can 
be analyzed at high precision. The evolution equations are ordinary differential equations, which 
can be solved analytically. The Wilson coefficients [5-8] and anomalous dimensions are given in 
terms of multiple nested harmonic sums [9,10], which obey algebraic relations [9-12]. Each finite 
harmonic sum can be represented in terms of a Mellin transform of a harmonic polylogarithm [13] 
weighted by l/(l±a;) and a polynomial of harmonic sums of lower weight. It turns out that for all 
basic functions at least up to weight five, as needed for the physical observables mentioned above, 
the harmonic polylogarithms are Nielsen integrals [14,15]. Furthermore, structural relations exist 
between various Mellin transforms of weighted Nielsen integrals [16]. They are due to relations 
at rational values of the Mellin variable, differentiation w.r.t. the Mellin variable, integration 
by parts relations and functional identities of Nielsen integrals. These relations, together with 
the quasi-shuffle algebra of harmonic sums allow to build the massless single-scale quantities 
in perturbative QED and QCD out of a low number of basic functions at a given order in 
perturbation theory. 

The anomalous dimensions and Wilson coefficients in Mellin space, for instance, can be 
written in terms of these basic functions and a number of their derivatives, in addition to 
polynomials including Euler’s ^-functions as well as rational functions of the Mellin variable for 
complex values of N. Then, on this basis, results in rr-space, are obtained by means of a single 
inverse Mellin transform, which has to be carried out numerically. For the class of functions 
contributing to the 2-loop Wilson coefficients the analytic continuations were given in Ref. [17]. 
Simple analytic continuations of nested sums are known for long, cf. e.g. Ref. [18]. This direction 
has been followed in Ref. [19] recently. However, in many cases the numerical convergence is 
very slow. 

To obtain an appropriate representation for the 3-loop anomalous dimensions and 2-loop 
Wilson coefficients, we exploit the algebraic and structural relations between the nested harmonic 
sums. For the 2-loop anomalous dimensions only one non-trivial function contributes aside from 
Euler’s ip^(N) functions [20,21], All known massless 2-loop Wilson coefficients can be expressed 
by four more functions [22,23]. For the 3-loop anomalous dimension the Mellin transforms 

M[/(x)](AT) = f dx x N f(x) (1.1) 

Jo 

of the following eight functions form the set of additional basic functions : 

LrQe) £1,3(3;) £2,2(3;) Li|(x) £2,2(~ff) ~ IA 2 2 {- x )/2 

1 ± a; ’ 1+x ’ l±x ’ 1 + x’ l±x ' [ } 

Here the functions Li n (a;) and S p . n {x) denote polylogarithms [24] and Nielsen-integrals [14,15], 
respectively. In terms of harmonic polylogarithms [13], they read 

Li 4 (a;) = -£0,0,0,1(37), £ 1 , 3 ( 3 :) = £0,1,1, 1 ( 3 ;), £ 2 , 2 ( 3 ;) = £0,0,1, 1 ( 3 ;), 

Li^) = (£o,i(3:)) 2 , S 2 , 2 (-x) = £ 0 , 0 ,- 1 ,- 1 ( 3 :), Lia(-3;) = (-H 0 ,-i(x)) 2 , (1.3) 

where the algebra can be applied to products of harmonic polylogarithms to express them in the 
basis of single harmonic polylogarithms. 
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In the present paper we derive fast and precise numerical representations of these Mellin 
transforms for complex values of N. Previously numerical representations of the 3-loop anoma¬ 
lous dimensions were given in [3,4], 1 Following earlier investigations [17], we aim at numerical 
representations for the individual basic functions, the building blocks of the massless Wilson 
coefficients and the anomalous dimensions, to high precision. Thus, in principle, the present 
approach is tunable to arbitrary accuracy. 


2 Mellin Transforms of the Type f(x)/( 1 + x ) 

As outlined in Ref. [17] the Mellin transforms containing the denominator 1/(1 + x) are given 
by: 


M 


' f( x ) ' 

1 + x 


(N) 


ln(2) • /(1) — f dx ln(l + x) [Nf(x ) + xf'(x)] 

Jo 

20 

ln(2) ■ /(l) - ]T “7 (JV M[/(i)](iV + k - 1) + M [/'(*)] (N + k)). 

k=1 


( 2 . 1 ) 


Accurate representations for the function ln(l + x) can be given using the minimax method for 
numerical approximation: 


20 

ln(l + x) ~ ^ aj^x fc 

k=1 


which we further improved to an accuracy of 3 x 10 16 . The coefficients read: 


A 

4 1 ’ 

4 1 ’ 

a^ 1} 

(i) 


a, 


9 

i (1) 

ni 

4.3 

i (1) 

l 15 

‘17 

l (1) 

‘l9 


0.9999999999999925D-0 a^ 
0.3333333332641123D-0 a ( { ] 
0.1999999561535526D-0 
0.1428517138099479D-0 a^ 
0.1109128496887138D-0 
0.8826572954250856D-1 a{^ 
0.6225829212455825D-1 
0.2890194939889559D-1 a{^ 
0.6003156359511387D-2 0$ 
0.3172112728405899D-3 ajj* 


—0.4999999999988568D-0 
—0.2499999977763199D-0 
—0.1666660875051348D-0 
—0.1249623936313475D-0 
—0.9918652787800788D-1 
—0.7643209265133132D-1 
—0.4572477090315515D-1 
—0.1496621145891488D-1 
—0.1731328252868496D-2 
—0.2760099875146713D-4 


For later use we also parameterize ln 2 (l + x) 


24 

ln 2 (l + x) ~ ^ , 

k=2 


delated earlier approximate representations were given in [25]. 


( 2 . 2 ) 


(2.3) 
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with 


(2) 

«2 

= 1.0000000000000000D-0 

(2) 

a 3 

—0.9999999999999985D-0 

(2) 

a\ 

= 0.9166666666663948D-0 

II 

-0.8333333333136118D-0 

(2) 

«6 

= 0.7611111103508889D-0 

4 2) = 

—0.6999999819735105D-0 

(2) 

«8 

= 0.6482139985629993D-0 

11 

<3 

—0.6039649964806160D-0 

(2) 

°10 

= 0.5657662306410356D-0 

2 

a n = 

—0.5323631718571445D-0 

(2) 

°12 

= 0.5024238774786239D-0 

(2) 

°13 = 

—0.4738508288315496D-0 

(2) 

°14 

= 0.4427472719775835D-0 

(2) 

a 15 = 

—0.4029142806330511D-0 

(2) 

°16 

= 0.3476841543351489D-0 

(2) 

a 17 = 

— 0.2748590021353420D-0 

a (2) 

a 18 

= 0.1915627642585285D-0 

a (2) - 

a 19 — 

-0.1130763066428224D-0 

(2) 

a 20 

= 0.5415661067306229D-1 

(2) 

«21 = 

— 0.1999877298940919D-1 

(2) 

a 22 

= 0.530362443938841ID-2 

(2) 

a 23 — 

— 0.8944156375768203D-3 

(2) 

«24 

= 0.7179502917974332D-4 




The following Mellin transforms contribute, cf. also Refs. [9,10,17], 


M [Li 4 (x)] (N - 1) 
M [Li^(x)] (N - 1) 
M[S , 1) 3(x)] (N- 1) 
M [S' h3 (x)] (IV-1) 
M [Li^x)] (IV-1) 

M [(Li|(*))l (IV - 1) 
M[S 2)2 (x)] (IV-1 ) 
(N- 1) 

M [S 2 , 2 (-x) - Ul(-x)/2] (IV-1) 


M 


(S 2 , 2 (-*)-Li 2 (-x)/2)'l (IV-1) 


C4 c 3 c 2 S 1 (N) 

N N 2 N 3 N 4 

M[Li 3 (x)] (IV-2) 

C4 (^i(iV)) 3 S 1 {N)S 2 {N) 

N 6 N 2 2N 2 

— ^ M [ln 3 (l — x)] (IV -2) 


S 3 (N) 

3N 2 


1 

N 


C2 2 _4Cs_ sw s*m 

s 'N s N N 


+—(( Sl (N)) 2 + S 2 (N)) 

—2 M [In(1 — x)Li 2 (x)] (At - 2) 

C4 
4 At 


C 3 (Si(N)) 2 + S 2 (N) 
N 2 + 2 At 3 


M[S 1)2 (x)] (At - 2) 


1 

N 
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(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 
(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


-C 2 2 + 2Li 4 (l/2) + -( 3 In(2) - -C 2 (hi(2)) 2 


+ -(ln(2)) 4 - M[S lt2 (-x)] (At — 1) 


1 

'TV 


r c 2 2 


+ M [In(1 + .r)l.i 2 (—.r)] (N - 1) 


M[5 1)2 (-x)] (At-2 ) 

+ M [In(1 + x)Li 2 (—x)] (At - 2) 


( 2 . 12 ) 

(2.13) 
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with 


M [Li 3 (x)] (N - 1) 
M [ln 3 (l - x)] (N - 1) 
M[S 1i2 (x)] (N- 1) 
M [S 1(2 (-x)] (N - 1) 


C 3 C 2 S 1 (N) 

N N 2 N 3 

(■ Si(N)) 3 0 S 1 (N)S 2 (N) S 3 (N ) 


N N 

C 3 QSi(iV)) 2 5 2 (iV) 
IV 2AT2 2AT2 

Ca (ln(2)) 2 (-1)” 

81V 2 N 2 N 2 


■2(M2)) 2 


IV 


(2.14) 

(2.15) 

(2.16) 


S' 1 ,_ 1 (7V)+ln(2) (S^N) - S-^N)) 


1 A a< 2 > 


M [In(1 — x)Li 2 (x)] (N — 1) = 


C 3 _ j_ \ ^ 

SAT 2N^'N + k 

k=2 

1 

N 


(2.17) 


-2Cs - C 2 S 1 W + ^ ((Si(JV)) 2 + S 2 (JV)) + S 2 ,,(l\r) 


(2.18) 


20 


M [ln(l + x)Li 2 (—x)] (IV — 1) ~ M[Li 2 (-:c)](lV-l + fc) 


(2.19) 


fc= 1 


M [Li 2 (—*)] (IV - 1) = -A + ^ - i-hl(5_i(Ar) + ln(2)) (2.20) 

The functions 5 'i j _i( 1V) and £ 2,1 (IV) have been parameterized in Ref. [17] 


S 2ll (N) = M 


where 


(-l^M 

M 

/ Li 2 (x) 


\ 1 — X 

M 

ln(l + x) 

1 + X 



ln(l + x) 


(■ N) ~ ln(2) (S 1 (N) - S-i(N)) + -(ln(2)) 2 (2.21) 


1 + x 

(A0 + c 2 Si (N), 


1 /V 24 u (2) 

(N)~-(in(2)) 2 --x; 


( 2 . 22 ) 


2 ^ N + k 

k=2 


(2.23) 


3 Mellin Transforms of the Type (/(a?)/( 1 — a?)) + 

Three functions of this type contribute. In case the numerator function is analytic at x — 1 one 
may represent the functions of this class using the minimax method directly. Both Li^x) and 
S 2j2 (x) have branch points at x — 1. Therefore we derive an analytic representation f t (x) for 
[fi(x) — fi( 1)] /(I — x) valid for the region x —» 1 and determine the remainder part using the 
minimax method: 


fj(x ) ~ /»(!) 

1 — x 


M%) + XI4°®* • 

fc=0 


(3.1) 
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For [Li 4 (x) — C 4 ] /(1 — x ) one obtains 


fi( x ) 


<3 + |(C 2 - 00(1 - x) + F ln(l - *) +1 -ii - | ] (! - x) 2 




a; 


— ln(l — x) + — (2 ~ — 
24 v ; 12 s 5 


599 


5 , \ 137 

16 ln(1 -" ) + 360 C2 


1440 
_79__Cs 

192 6 


(1 

(1 


x) 


x) 


(3.2) 


c 

c 

c 

c 

c 

c 


( 1 ) 

0 

( 1 ) 

2 

( 1 ) 

4 

( 1 ) 

6 

( 1 ) 

8 

( 1 ) 

10 


c 


( 1 ) 

12 


—3.187045493829754D-1 cj 1} 
—3.926780960319761D-0 
—2.764070067739643D-0 
—1.412597571664758D-2 
2.072912823276118D-1 
3.100508803799690D-1 cj 1 ^ 
7.682650942255444D-2 


1.752102582962004D-0 
4.533622455411171D-0 
7.705635337822153D-1 

— 1.024735147728843D-1 
—3.094485142894180D T 
—2.017510797419543D-1 

— 1.310258217741916D-2 


with k mSLX = 13. The Mcllin transform of fi{x) is given by a linear combination of 


M [(1 - x) M ] (N) = B(N+1,M+1) (3.3) 

M [(1 -x) M ln(l -x)] (N) = [ip(M + 1) — i(j(M + N + 2)\ B(N + l, M + 1) (3.4) 

M [(1 -x) M ln 2 (l -x)j (N) = + 1) — + IV + 2)] 

+ [rj){M + 1) - ip(M + N + 2)} 2 }B(N + 1, M + 1) . (3.5) 


h ( x ) = (1 — x ) (^-^ 5 + 


/ 947 
’ \3600 


x 


720" 

11683 4 
- x 4 + 


ouu 


15221 

- x 


4827 2 
-x 2 + 


18511 


+ 


N 7200 3600 800 3600 

104641 | 11675141 15334867^ | 14820287 

504000^ + 9072000 X ~ 4536000 X + 3024000 
19680697 29645831 

4536000 X + 1134000 
43 , 667 


- x - x -1 - x -I mil 

3600 800 3600 150 J v 


7 


x 6 + -X- 


- X - Jb - - Jb - Jb - - 

210 420 140 140 140 


2341 .. 
-x 6 


853 2 

-x 


+ 


617 363 N 


-x 


C3 


5 



II 

So 

0 

2.319102959447070D-1 

c (2) - 
c i — 

— 1.341835247346528D-0 

Cb 

II 

3.141213505948200D-0 

c (2) _ 
c 3 

—3.689298402805891D-0 

C (2) - 
c 4 — 

2.066708069184852D-0 

c (2) - 

c 5 

—0.221726622455410D-0 

C (2) - 

C 6 

—2.866874606436715D-1 

0 

II 

1.886633904893633D-1 

c (2) - 

c 8 

—5.104536639955758D-1 

0 

03 s 

II 

1.814194435712243D-0 

J 2) _ 
c 10 — 

—4.810364434281610D-0 

C (2) - 

C 11 — 

9.670512224938253D-0 

C (2) - 

c 12 ~ 

— 14.74512233105934D-0 

c (2) - 

c 13 — 

16.93126634317140D-0 

c (2) - 

c 14 — 

— 14.39776586765419D-0 

c (2) - 

8.789666564310812D-0 

c (2) - 

c 16 — 

—3.642212170527524D-0 

c (2) - 

c 17 — 

9.173661876074224D-1 

c (2) ~ 

l 18 — 

— 1.060348165292345D-1 
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The Nielsen integrals S np {x) are analytic at x = —1. Therefore [*S' 2 > 2 (— x) — Li|(—a ?)/2 — ci] / 
(1 — x) can be represented by a polynomial directly, where 

c i = _ g^2 2 + 2Li 4 (l/2) + -( 3 ln(2) - -C 2 (ln(2 )) 2 + —(ln(2 )) 4 . (3.7) 

We will do this for the functions [S 2j 2 {~x) — c 4 — C 2 2 /8] /(1 — x) and [Li^ — x) — C 2 2 / 4 ] /(I — x) 
separately. 

The following representation for M[( 5 ' 2 ) 2 (— x) — c 4 — C 2 2 / 8 )/(1 — x)] is obtained: 


M 


S2,2(—%) ~ Cl — C2 2 /8 

1 — x 


(iV) 


E 



IV + k + 1 


5 


6q 1} = —0.8778567156865530D T b[ l) 

b ( ;p = 0.3721432843134470D-1 

b ( { ] = 0.1030460620911406D-1 lx t } 

b ( b L) = 0.4208927186604172D-2 ll l) 

6g 1} = 0.2130242091464232D-2 

6^ = 0.1230894254048324D-2 b^ 

b[l] = 0.7755737888131504D-3 b\'] 

b[ l J = 0.5031055861991529D-3 

b[\} = 0.2868111320615437D-3 b$ 

b{^ = 0.1068980760727356D-3 6$ 

6 2 q = 0.1797845625225957D-4 b^ 

b$ = 0.7830613264154134D-6 b$ 


—0.8778567156865530D-1 
—0.1834122712421077D-1 
—0.6362060457158173D-2 
—0.2933929772309200D-2 
—0.1597936899156854D-2 
—0.9689535635431825D-3 
—0.6263517402742744D-3 
—0.3922594818742940D-3 
—0.1887234737089442D-3 
—0.4971730708906830D-4 
—0.4695889540721375D-5 
—0.6232207394074941D-7 


The corresponding polynomial in x has an accuracy of 2 x 10 15 in [0,1]. 
Similarly, for M[(Li 2 ( —x) — C 2 2 /4)/(l — x)](N) one obtains: 


M 


Li 2 


-x 


C2 2 /4 


1 — x 


(N) = 


11 


1 ,( 2 ) 


^ N + k + l 

k =0 


(3.8) 


(3.9) 


6 



with 


(2) 

0 — 

—0.6764520210934552D-0 

6< 2) = 

— 0.6764520137562308D-0 

(2) _ 

2 

0.3235476094265664D-0 

<$> = 

—0.1764446743143206D-0 

(2) 

4 — 

0.1081940672246993D-0 

4 2 > = 

—0.7181309059958118D-1 

(2) 

6 “ 

0.4940999469881481D-1 

4 2 ’ = 

—0.3290941711692155D-1 

(2) 

8 

0.1916664887064280D-1 

6 (2 > - 
°9 

—0.8589741767655388D-2 

(2) 

10 — 

0.2508898780543465D-2 

b {2) ~ 

°11 ~ 

— 0.3476710199486832D-3 


The corresponding polynomial in x has an accuracy of 2.5 x 10 11 for x e [0.1]. 

In Table 1 the values of the maximum relative errors are given for the first 40 Mcllin moments. 
Depending on the function accuracies better than 9 • 10~ 10 to 8 • 10” 13 are obtained. We also 
compared the representation of the basic functions through numerical inversion of the Mcllin 
moments by a complex contour integral. In the range x e [10 -6 , 0.98] maximal relative errors of 
2 to 8 • 10 -8 were obtained, irrespective of the value of the basic functions. In some regions in 
x even up to several order of magnitude better results were obtained. 



Moments, IV < 40 

Inversion 

[Li 4 (x) — C 4 ]/(l — a?) 

4.0D-12 

5.1D-8 

[£ 2 , 2 (t) — C 4 / 4]/(1 — x) 

8.2D-10 

5.4D-8 

[^(-z) - Ci]/(1 - x) 

1.8D-10 

4.5D-8 

m-x)-Q/A}/(i-x) 

3.0D-10 

4.3D-8 

Li 4 (x)/(1 + x ) 

1.2D-12 

8.0D-8 

S\^{x)/ (1 + x) 

2.2D-12 

3.3D-8 

S 2 fl{x)/{l+x) 

1.1D-12 

1.8D-8 

Li 2 (»/(! + u) 

2.2D-12 

3.0D-8 

[S 2 , 2 (-x) - Lio(—x)/2 ]/(l + x) 

7.2D-13 

3.2D-8 


Table 1: Relative accuracies of the representations. 


4 Conclusions 

Single scale quantities in massless QED and QCD such as the Wilson coefficients and anomalous 
dimensions can be represented by nested harmonic sums. Their representation may be reduced 
to a small set of basic functions in addition to Euler’s i/j( k \N) functions, each of the basic func¬ 
tions corresponding to one particular single harmonic sums. This is possible due to algebraic and 
structural relations between the harmonic sums. Up to the level of the 3-loop anomalous dimen¬ 
sions 13 such basic functions occur. We have derived fast and precise numerical representations 
for the eight new functions contributing at the level of the 3-loop anomalous dimensions in QCD. 
In the region x e [10~ 6 , 0.98] which is physically most interesting their relative accuracy is better 
than 10 -7 . If needed, the procedure outlined in the present paper can readily be generalized to 
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a higher level of accuracy at the expense of more coefficients in the representation. Due to use of 
polynomial forms in the analytic continuations very fast representations are obtained well suited 
for the use in QCD evolution codes. The basic functions discussed in the present paper form 
a part of the building blocks for the 3-loop coefficient functions and other higher order single 
scale quantities. 

FORTRAN routines are available from the authors upon request. 

Acknowledgment. 

This paper was supported in part by DFG Sonderforschungsbereich Transregio 9, Comput- 
ergestiitzte Theoretische Physik. S.M. acknowledges the hospitality of the FNAL Theory Divi¬ 
sion during completion of the work. 


References 

[1] M. Gluck, E. Reya and A. Vogt, Z. Phys. C48 (1990) 471. 

[2] J. Bliimlein and A. Vogt, Phys. Rev. D58 (1998) 014020, hep-ph/9712546. 

[3] S. Moch, J.A.M. Vermaseren and A. Vogt, Nucl. Phys. B688 (2004) 101, hep-ph/0403192. 

[4] A. Vogt, S. Moch and J.A.M. Vermaseren, Nucl. Phys. B691 (2004) 129, hep-ph/0404111. 

[5] D. 1. Kazakov and A. V. Kotikov, Nucl. Phys. B307 (1988) 721 [Erratum-ibid. B345 (1990) 
299], 

[6] J. Sanchez Guillen, J. Miramontes, M. Miramontes, G. Parente and O. A. Sampayo, Nucl. 
Phys. B353 (1991) 337. 

[7] W.L. van Neerven and E.B. Zijlstra, Phys. Lett. B272 (1991) 127; B273 (1991) 476; B297 
(1992) 377; Nucl. Phys. B383 (1992) 525. 

[8] S. Moch and J.A.M. Vermaseren, Nucl. Phys. B573 (2000) 853, hep-ph/9912355. 

[9] J.A.M. Vermaseren, Int. J. Mod. Phys. A14 (1999) 2037, hep-ph/9806280. 

[10] J. Bliimlein and S. Kurth, Phys. Rev. D60 (1999) 014018, hep-ph/9810241. 

[11] J.M. Borwein and R. Girgensohn, Electron. J. Combinatorics 3 (1996) R23 (Appendix by 
D.J. Broadhnrst), CECM research report 95-053. 

[12] J. Bliimlein, Comput. Phys. Commun. 159 (2004) 19, hep-ph/0311046. 

[13] E. Remiddi and J.A.M. Vermaseren, Int. J. Mod. Phys. A15 (2000) 725, hep-ph/9905237. 

[14] N. Nielsen, Nova Acta Leopoldina (Halle) 90 (1909) 123. 

[15] K.S. Kolbig, SIAM J. Math. Anal. 17 (1986) 1232. 

[16] J. Bliimlein, Nucl. Phys. Proc. Suppl. 135 (2004) 225, hep-ph/0407044. 

[17] J. Bliimlein, Comput. Phys. Commun. 133 (2000) 76, hep-ph/0003100. 

[18] J. Bliimlein and S. Kurth, (1997), hep-ph/9708388. 



[19] A.V. Kotikov and V.N. Velizhanin, (2005), hep-ph/0501274. 

[20] A. Gonzalez-Arroyo, C. Lopez and F.J. Yndnrain, Nucl. Phys. B153 (1979) 161; 

[21] A. Gonzalez-Arroyo and C. Lopez, Nucl. Phys. B166 (1980) 429. 

[22] J. Bliimlein and S. Moch, to appear. 

[23] J. Bliimlein and V. Ravindran, (2005), hep-ph/0501178. 

[24] L. Lewin, Polylogarithms and Associated Functions (North Holland, New York, 1981), ISBN 
0-444-00550-1. 

[25] W.L. van Neerven and A. Vogt, Nucl. Phys. B568 (2000) 263, hep-ph/9907472; B588 
(2000) 345, hep-ph/0006154; B603 (2001) 42, hep-ph/0103123; Phys. Lett. B490 (2000) 
111, hep-ph/0007362. 


9 



